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Abstract 

Physical systems exhibiting stochastic or chaotic behavior are often amenable 
to treatment by random matrix models. In deciding on a good choice of model, 
, random matrix physics is constrained and guided by symmetry considerations. 

\ The notion of 'symmetry class' (not to be confused with 'universality class') 

expresses the relevance of symmetries as an organizational principle. Dyson, 
in his 1962 paper referred to as The Threefold Way, gave the prime classifica- 
tion of random matrix ensembles based on a quantum mechanical setting with 
symmetries. In this article we review Dyson's Threefold Way from a modern 
perspective. We then describe a minimal extension of Dyson's setting to in- 
corporate the physics of chiral Dirac fermions and disordered superconductors. 
In this minimally extended setting, where Hilbert space is replaced by Fock 
space equipped with the anti- unitary operation of particle- hole conjugation, 
symmetry classes are in one-to-one correspondence with the large families of 
Riemannian symmetric spaces. 



1 Introduction 

In Chapter 2 of this handbook^], the historical narrative by Bohigas and Wei- 
denmiiller describes how random matrix models emerged from quantum physics, 
more precisely from a statistical approach to the strongly interacting many- 
body system of the atomic nucleus. Although random matrix theory is nowa- 
days understood to be of relevance to numerous areas of physics, mathematics, 
and beyond, quantum mechanics is still where many of its applications lie. 

1 The present article is to be Chapter 3 of the Oxford Handbook of Random Matrix Theory. 
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Quantum mechanics also provides a natural framework in which to classify 
random matrix ensembles. 

In this thrust of development, a symmetry classification of random matrix 
ensembles was put forth by Dyson in his 1962 paper The Threefold Way: alge- 
braic structure of symmetry groups and ensembles in quantum mechanics, where 
he proved (quote from the abstract of |Dys62|) "that the most general matrix 
ensemble, defined with a symmetry group which may be completely arbitrary, 
reduces to a direct product of independent irreducible ensembles each of which 
belongs to one of the three known types". The three types known to Dyson 
were ensembles of matrices which are either complex Hermitian, or real symmet- 
ric, or quaternion self-dual. It is widely acknowledged that Dyson's Threefold 
Way has become fundamental to various areas of theoretical physics, includ- 
ing the statistical theory of complex many-body systems, mesoscopic physics, 
disordered electron systems, and the field of quantum chaos. 

Over the last decade, a number of random matrix ensembles beyond Dyson's 
classification have come to the fore in physics and mathematics. On the physics 
side these emerged from work [Ver94| on the low-energy Dirac spectrum of 
quantum chromodynamics, and also from the mesoscopic physics of low-energy 
quasi-particles in disordered superconductors [AZ97J. In the mathematical re- 
search area of number theory, the study of statistical correlations of the values 
of Riemann zeta and related L-functions has prompted some of the same gen- 
eralizations |KS99j. It was observed early on |AZ97j that these post-Dyson 
ensembles, or rather the underlying symmetry classes, are in one-to-one corre- 
spondence with the large families of symmetric spaces. 

The prime emphasis of the present handbook article will be on describing 
Dyson's fundamental result from a modern perspective. A second goal will be to 
introduce the post-Dyson ensembles. While there seems to exist no unanimous 
view on how these fit into a systematic picture, here we will follow [HHZ05J to 
demonstrate that they emerge from Dyson's setting upon replacing the plain 
structure of Hilbert space by the more refined structure of Fock spaceJl The 
reader is advised that some aspects of this story are treated in a more leisurely 
manner in the author's encyclopedia article |Zir04| . 

To preclude any misunderstanding, let us issue a clarification of language 
right here: 'symmetry class' must not be confused with 'universality class' ! 
Indeed, inside a symmetry class as understood in this article various types of 
physical behavior are possible. (For example, random matrix models for weakly 
disordered time-reversal invariant metals belong to the so-called Wigner-Dyson 
symmetry class of real symmetric matrices, and so do Anderson tight-binding 
models with real hopping and strong disorder. The former are believed to 

2 We mention in passing that a classification of Dirac Hamiltonians in two dimensions has 
been proposed in [BL02[ . Unlike ours, this is not a symmetry classification in Dyson's sense. 
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exhibit the universal energy level statistics given by the Gaussian Orthogo- 
nal Ensemble, whereas the latter have localized eigenfunctions and hence level 
statistics which is expected to approach the Poisson limit when the system size 
goes to infinity.) For this reason the present article must refrain from writing 
down explicit formulas for joint eigenvalue distributions, which are available 
only in certain universal limits. 

2 Dyson's Threefold Way 

Dyson's classification is formulated in a general and simple mathematical set- 
ting which we now describe. First of all, the framework of quantum theory 
calls for the basic structure of a complex vector space V carrying a Hermitian 
scalar product (•,•): V x V — > C (Dyson actually argues |Dys62| in favor of 
working over the real numbers, but we will not follow suit in this respect.) For 
technical simplicity, we do join Dyson in taking V to be finite-dimensional. In 
applications, V — C n will usually be the truncated Hilbert space of a family of 
disordered or quantum chaotic Hamiltonian systems. 

The Hermitian structure of the vector space V determines a group U(V) 
of unitary transformations of V. Let us recall that the elements g £ U(V) are 
C-linear operators satisfying the condition (gv,gv') = (v,v ! ) for all v,v' £ V. 

Building on the Hermitian vector space V, Dyson's setting stipulates that 
V be equipped with a unitary group action 

G xV^V, (g,v) ^ Pv (g)v , p v (g)eTJ(V). (2.1) 

In other words, there is some group Go whose elements g are represented on V 
by unitary operators pvig)- This group Go is meant to be the group of joint 
(unitary) symmetries of a family of quantum mechanical Hamiltonian systems 
with Hilbert space V. We will write pv(g) — g f° r short. 

Now, not every symmetry of a quantum system is of the canonical unitary 
kind. The prime counterexample is the operation, T, of inverting the time 
direction, called time reversal for short. It is represented on Hilbert space V by 
an anti-unitary operator T = py(T), which is to say that T is complex anti- 
linear and preserves the Hermitian scalar product up to complex conjugation: 

T(zv) = zTv, (Tv,Tv') = (v,v>) (z G C ; v, v £ V) . (2.2) 

Another operation of this kind is charge conjugation in relativistic theories such 
as the Dirac equation for the electron and its anti-particle, the positron. 

Thus in Dyson's general setting one has a so-called symmetry group G = 
Go U Gi where the subgroup Go is represented on V by unitaries, while G\ (not 
a group) is represented by anti-unitaries. By the definition of what is meant 
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by a 'symmetry', the generator of time evolution, the Hamiltonian H, of the 
quantum system is fixed by conjugation gHg^ 1 = H with any g E G. 

The set G\ may be empty. When it is not, the composition of any two 
elements of G\ is unitary, so every g E G\ can be obtained from a fixed element 
of G\ , say T, by right multiplication with some U E Go : g = TU. The 
same goes for left multiplication, i.e., for every g E G\ there also exists V E 
Go so that g = U'T. In other words, when G± is non-empty, Go C G is a 
proper normal subgroup and the factor group G/Gq ~ Z2 consists of exactly 
two elements, Go and TGo = G\ . For future use we record that conjugation 
U I— > TUT^ 1 =: a(U) by time reversal is an automorphism of Go • 

Following Dyson |Dys62| we assume that the special element T represents 
an inversion symmetry such as time reversal or charge conjugation. T must 
then be a (projective) involution, i.e., T 2 = z x Idy with 7^ z E C , so that 
conjugation by T 2 is the identity operation. Since T is anti- unitary, z must 
have modulus \z\ = 1 , and by the C-antilinearity of T the associative law 

zT = T 2 ■ T = T ■ T 2 = Tz = zT (2.3) 

forces z to be real, which leaves only two possibilities: T 2 = ±Idy . 

Let us record here a concrete example of some historical importance: the 
Hilbert space V might be the space of totally anti-symmetric wave functions 
of n particles distributed over the shell-model space of an atom or an atomic 
nucleus, and the symmetry group G might be G = O3 U TO3 , the full rotation 
group O3 (including parity) together with its translate by time reversal T. 

In summary, Dyson's setting assumes two pieces of data: 

• a finite-dimensional complex vector space V with Hermitian structure, 

• a group G = Go U TGo acting on V by unitary and anti- unitary operators. 

It should be stressed that, in principle, the primary object is the Hamiltonian, 
and the symmetries G are secondary objects derived from it. However, adopting 
Dyson's standpoint we now turn tables to view the symmetries as fundamental 
and given and the Hamiltonians as derived objects. Thus, fixing any pair (V, G) 
our goal is to elucidate the structure of the space of all compatible Hamiltoni- 
ans, i.e., the self-adjoint operators H on V which commute with the G-action. 
Such a space is reducible in general: the G-compatible Hamiltonian matrices 
decompose as a direct sum of blocks. The goal of classification is to enumerate 
the irreducible blocks that occur in this setting. 

While the main objects to classify are the spaces of compatible Hamiltonians 
H, we find it technically convenient to perform some of the discussion at the 
integrated level of time evolutions Ut = Q~ itH / h instead. This change of focus 
results in no loss, as the Hamiltonians can always be retrieved by linearization 
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in t at t = . The compatibility conditions for U = Ut read 

U = g Ug Q 1 =g l U- 1 g^ 1 (for all g a G G a ) . (2.4) 

The strategy will be to make a reduction to the case of the trivial group 
Go = {Id}. The situation with trivial Go can then be handled by enumeration 
of a finite number of possibilities. 

2.1 Reduction to the case of G = {Id} 

To motivate the technical reduction procedure below, we begin by elaborating 
the example of the rotation group O3 acting on a Hilbert space of shell-model 
states. Any Hamiltonian which commutes with Go = O3 conserves total an- 
gular momentum, L, and parity, 7T, which means that all Hamiltonian matrix 
elements connecting states in sectors of different quantum numbers (L, it) van- 
ish identically. Thus, the matrix of the Hamiltonian with respect to a basis of 
states with definite values of (L, it) has diagonal block structure. 03-symmetry 
further implies that the Hamiltonian matrix is diagonal with respect to the 
orthogonal projection, M, of total angular momentum on some axis in position 
space. Moreover, for a suitable choice of basis the matrix will be the same for 
each M- value of a given sector (L,tt). 

To put these words into formulas, we employ the mathematical notions of 
orthogonal sum and tensor product to decompose the shell-model space as 

v - *W)> V (L)7r) =C^)®C 2i+1 , (2.5) 

L>0; 7T=±1 

where m(L, it) is the multiplicity in V of the 03-representation with quantum 
numbers (L, it). The statement above is that all symmetry operators and com- 
patible Hamiltonians are diagonal with respect to this direct sum, and within 
a fixed block V(£ j7r ) the Hamiltonians act on the first factor C m( - L '^ and are 
trivial on the second factor C 2L+1 of the tensor product, while the symmetry 
operators act on the second factor and are trivial on the first factor. Thus 
we may picture each sector V(l,-k) as a rectangular array of states where the 
Hamiltonians act, say, horizontally and are the same in each row of the array, 
while the symmetries act vertically and are the same in each column. 

This concludes our example, and we now move on to the general case of 
any group Go acting reductively on V. To handle it, we need some language 
and notation as follows. A Go-representation X is a C-vector space carrying a 
Go-action Go x X — > X by (g, x) i->- px(g)x. If X and Y are Go-representations, 
then by the space HomQ, (X, Y) of Go-equivariant homomorphisms from X to 
Y one means the complex vector space of C-linear maps ip : X — > Y with 
the intertwining property py^g)^ = tpPx(g) f° r all g G Go . If X is an ir- 
reducible Go-representation, then Schur's lemma says that Homc (X, X) is 
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one-dimensional, being spanned by the identity, Idx • For two irreducible Go- 
representations X and Y, the dimension of Hom<3 (X, Y) is either zero or one, 
by an easy corollary of Schur's lemma. In the latter case X and Y are said to 
belong to the same isomorphism class. 

Using the symbol A to denote the isomorphism classes of irreducible Go- 
representations, we fix for each A a standard representation space R\ . Note that 
dim HoniG (R\ , V ) counts the multiplicity in V of the irreducible representation 
of isomorphism class A . In our shell-model example with Go = O3 we have 
A= (L,vr), R X = C 2L+1 , and dim Homg (R\ , V) =m(L,vr). 

The following statement can be interpreted as saying that the example ad- 
equately reflects the general situation. 

Lemma 2.1 Let Go act reductively on V . Then 



is a Go-equivariant isomorphism. 

Remark. The decomposition offered by this lemma perfectly separates the 
unitary symmetry multiplets from the dynamical degrees of freedom and thus 
gives an immediate view of the structure of the space of Go-compatible Hamil- 
tonians. Indeed, the direct sum over isomorphism classes (or 'sectors') A is 
preserved by the symmetries Go as well as the compatible Hamiltonians H; 
and Go is trivial on Homg (R\ , V) while the Hamiltonians are trivial on R\ . 

Next, we remove the time-evolution trivial factors R\ from the picture. To 
do so, we need to go through the step of transferring all given structure to the 
spaces E\ := rIomG (i?^ , V). 

2.1.1 Transfer of structure. 

We first transfer the Hermitian structure of V. In the present setting of a 
unitary Go-action, the Hermitian scalar product of V reduces to a Hermitian 
scalar product on each sector of the direct-sum decomposition of Lemma 12.11 
by orthogonality of the sum. Hence, we may focus attention on a definite sector 
E <g) R = E\ <g> R\ . Fixing a Go-invariant Hermitian scalar product (-,-)r on 
R = R\ we define such a product {-,-)e '■ E x E — > C by 



which is easily checked to be independent of the choice of r G R , r 7^ . 

Before carrying on, we note that for any Hermitian vector space V there 
exists a canonically defined C-antilinear bijection Cy : V — > V* to the dual 
vector space V* by Cy(v) := (v, -}y . (In Dirac's language this is the conversion 
from 'kef vector to 'bra' vector.) By naturalness of the transfer of Hermitian 
structure we have the relation Ce®r = Ce® Cr . 



nomG (R x ,V)®Rx^V, 



(V>A ® r\) ' y V tp\{r\) 



(tp,ip')E ■= {ip(r),if)'(r)) v /(r,r) R 



(2.6) 
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Turning to the more involved step of transferring time reversal T, we begin 
with a preparation. If L : V — > W is a linear mapping between vector spaces, we 
denote by L* : W* — > V* the canonical transpose defined by (L t f)(v) = f(Lv). 
Let now V be our Hilbert space with ket-bra bijection C = Cy . Then for any 
g £ U(V) we have the relation CgC~ 1 = (g~ 1 ) t because 

C(gv) = (gv,-) = (v^g- 1 -) = {g-^Civ) (v e V) . (2.7) 

Moreover, recalling the automorphism Go B g i— > a(g) = TgT" 1 of Go we obtain 

CT g = aig-y CT (g € G ) . (2.8) 

Thus, since C and T are bijective, the C-linear mapping CT : V — > V* is a 
Go-equivariant isomorphism interchanging the given Go-representation on V 
with the representation on V* by g H ?■ a{g^ 1 ) 1 . In particular, it follows that T 
stabilizes the decomposition V = ©a^A — ©a(-E'a ® R\) of Lemma |2~T1 

If T exchanges different sectors Va , the situation is very easy to handle (see 
below). The more challenging case is TV\ = V\ , which we now assume. 

Lemma 2.2 Let TV\ = V\ . Under the isomorphism V\ — E\ © R\ the time- 
reversal operator transfers to a pure tensor 

T = a © , a: E x ^ E x , : R x -»■ R X , 

with anti-unitary a and (3 . 

Proof. Writing E\ = E and Rx = R for short, we consider the transferred 
mapping CT : E © R — >■ £/* © ii* , which expands as GT = ]P <j){ © ^ with C- 
linear mappings 4>i . Since GT is known to be a Go-equivariant isomorphism, 
so is every map ipi : i? — > R*. By the irreducibility of R and Schur's lemma, 
there exists only one such map (up to scalar multiples). Hence CT is a pure 
tensor: CT = cp <g> ?/> . Using G = Ce®r = Ce © Gr we obtain T = a <S> (3 
with C-antilinear a = C~^ l cj) and /3 = C^xfj . Since the tensor product lets 
you move scalars between factors, the maps a and f3 are not uniquely defined. 
We may use this freedom to make f3 anti-unitary. Because T is anti-unitary, it 
then follows from the definition (j2.6p of the Hermitian structure of E that a is 
ant i- unitary. 

Remark. By an elementary argument, which was spelled out for the anti- 
unitary operator T in Eq. (|2.3|) . it follows that a 2 = eoTd^ and f3 2 = e^Id/j 
with e a , ep G {±1}- Writing T 2 = eyldy we have the relation e a e^ = ex ■ Thus 
when ep = — 1 the parity e a = — £t of the transferred time- reversal operator a 
is opposite to that of the original time reversal T. 

This change of parity occurs, e.g., in the case of Go = SU2 . Indeed, let 
R = R n be the irreducible SU2-representation of dimension n + 1 . It is a 
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standard fact of representation theory that R n is SU2-equivariantly isomorphic 
to i?* by a symmetric isomorphism ip = tp t for even n and skew-symmetric 
isomorphism tp = —i(i t for odd n . From (—l) n ip t = ip = Cr j3 and 

iP(v){v') = (/3v, v') R = W 2 v,(3v') R = </V, p 2 v) R = i>{v'){p 2 v) , (2.9) 

we conclude that /3 2 = (— l) n Id^ n . 

2.2 Classification 

By the decomposition of Lemma l2. li the space Zu(y)(Go) of Go-compatible time 
evolutions in U(V) is a direct product of unitary groups, 

z u(v) (G )^n A u ^- (2 - io) 

We now fix a sector V\ — E\ ® R\ and run through the different situations 
(of which there exist three, essentially) due to the absence or presence of a 
transferred time-reversal symmetry a : E\ — > E\ . 

2.2.1 Class A 

The first type of situation occurs when the set G\ of anti-unitary symmetries is 
either empty or else maps V\ ~ E\(&R\ to a different sector V\> , X^= X'. In both 
cases, the G-compatible time-evolution operators restricted to V\ constitute a 
unitary group U(E\) ~ Uat with N = dimE\ being the multiplicity of the 
irreducible representation R\ in V. The unitary groups \Jn or to be precise, 
their simple parts SU^ , are symmetric spaces (cf. Section 13. 4p of the A family 
or A series in Cartan's notation - hence the name Class A. In random matrix 
theory, the Lie group LV equipped with Haar measure is commonly referred to 
as the Circular Unitary Ensemble, CUE^r |Dys62a| . 

The Hamiltonians H in Class A are represented by complex Hermitian NxN 
matrices. By putting a LV-invariant Gaussian probability measure 

N 

d^{H) = c e- TrH2/2eT2 dH , dH = JJ dH u JJ dH jk dH kj , (2.11) 

i=l j<k 

on that space, one gets what is called the GUE - the Gaussian Unitary Ensem- 
ble - defining the Wigner-Dyson universality class of unitary symmetry. An 
important physical realization of that class is by electrons in a disordered metal 
with time-reversal symmetry broken by the presence of a magnetic field. 

2.2.2 Classes Al and All 

We now turn to the cases where T is present and TV\ = V\ ~ E\ ® R\ . We 
abbreviate E\ = E. From Lemma 12. II we know that T = a ® (3 is a pure tensor 
with anti- unitary a , and we have a 2 = e a Ue with parity e a = ex ■ 
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Using conjugation by a. to define an automorphism 

r : U(E) ->■ \J(E) , Ji^am" 1 , (2.12) 

we transfer the conditions (|2,4|) to V\ and describe the set Z-^ie)[G) of G- 
compatible time evolutions in U(-E) as 

Zv(E)(G) = {xe V(E) | r(x) = x- 1 } . (2.13) 

Now let U = U(-E) for short and denote by K C U the subgroup of r-fixed 
elements k = r(k) G 17. The set Zjj(G) is analytically diffeomorphic to the 
coset space U /K by the mapping 

U/K Zt/(G) Cf/, uif i-> mt^" 1 ) , (2.14) 

which is called the Cartan embedding of U/K into U. The remaining task is 
to determine K. This is done as follows. 

Recalling the definition Ceol = 4> an d using (Tg/c = (k~ 1 ) t CE we express 
the fixed-point condition k = r(k) = aka~ l as {k~ 1 ) 1 = (frkcf) -1 or, equivalently, 
(f) = k l (j)k , which means that the bilinear form associated with <p : E — >■ -E* , 

: E x E ^ C , (e, e') (-»• 0(e)(e') , (2.15) 

is preserved by fc G -K". By running the argument around Eq. (|2.9[) in reverse 
order (with the obvious substitutions f/? — )■ and /3 — >• a), we see that the 
non-degenerate form is symmetric if e a = +1 and skew if e a = — 1 . In the 
former case it follows that K = 0(E) ~ Oat is an orthogonal group, while in 
the latter case, which occurs only if N G 2N , K = USp(E) ~ USp^y is unitary 
symplectic. In both cases the coset space U/K is a symmetric space (cf. Section 
13. 4|) - a fact first noticed by Dyson in |Dys70| . 

Thus in the present setting of TV\ = V\ we have the following dichotomy 
for the sets of G-compatible time evolutions Zu^(G) ~ U/K : 

Class AI: U/K~\1 N /0 N {e a = +1) , 

Class All: U/K ~ Ujv/US PiV (e a = -1 , N G 2N) . 

Again we are referring to symmetric spaces by the names they - or rather their 
simple parts SUn/SOn and SUat/USp^ - have in the Cartan classification. 
In random matrix theory, the symmetric space Uat/Ojv (or its Cartan embed- 
ding into Uat as the symmetric unitary matrices) equipped with UAr-invariant 
probability measure is called the Circular Orthogonal Ensemble, COEjy , while 
the Cartan embedding of Uat/USp^ equipped with UAr-invariant probability 
measure is known as the Circular Symplectic Ensemble, CSE^v [Dys62a| . (Note 
the confusing fact that the naming goes by the subgroup which is divided out.) 
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Examples for Class AI are provided by time-reversal invariant systems with 
symmetry Go = (SU2) sp in • Indeed, by the fundamental laws of quantum physics 
time reversal T squares to (— 1) 2S times the identity on states with spin S. Such 
states transform according to the irreducible SU2-representation of dimension 
2S + 1, and from /3 2 = (— 1) 2S (see the Remark after Lemma 12. ip it follows 
that e a = = (— 1) 2S (— 1) = +1 in all cases. A historically important 

realization of Class AI is furnished by the highly excited states of atomic nuclei 
as observed by neutron scattering just above the neutron threshold. 

By breaking SU2-symmetry (i.e., by taking Go = {Id}) while maintaining 
T-symmetry for states with half-integer spin (say single electrons, which carry 
spin S = 1/2), one gets e a = ex = (— l) 25 = —1, thereby realizing Class All. An 
experimental observation of this class and its characteristic wave interference 
phenomena was first reported in the early 1980's [Bcr84] for disordered metallic 
magnesium films with strong spin-orbit scattering caused by gold impurities. 

The Hamiltonians H, obtained by passing to the tangent space of U/K at 
unity, are represented by Hermitian matrices with entries that are real num- 
bers (Class AI) or real quaternions (Class ^411). The simplest random matrix 
models result from putting iT-invariant Gaussian probability measures on these 
spaces; they are called the Gaussian Orthogonal Ensemble and Gaussian Sym- 
plectic Ensemble, respectively. Their properties delineate the Wigner-Dyson 
universality classes of orthogonal and symplectic symmetry. 

3 Symmetry Classes of Disordered Fermions 

While Dyson's Threefold Way is fundamental and complete in its general Hilbert 
space setting, the early 1990's witnessed the discovery of various new types of 
strong universality, which were begging for an extended scheme: 

• The introduction of QCD-motivated chiral random matrix ensembles (re- 
viewed by Verbaarschot in Chapter 32 of this handbook) mimicked Dyson's 
scheme but also transcended it. 

• Number theorists had introduced and studied ensembles of L-functions 
akin to the Riemann zeta function (see the review by Keating and Snaith 
in Chapter 24 of this handbook). These display random matrix phenom- 
ena which are absent in the classes A, AI, or All. 

• The proximity effect due to Andreev reflection, a particle-hole conversion 
process in mesoscopic hybrid systems involving metallic as well as su- 
perconducting components, was found |AZ97] to give rise to post-Dyson 
mechanisms of quantum interference (cf. Chapter 35 by Beenakker). 

By the middle of the 1990's, it had become clear that there exists a unifying 
mathematical principle governing these post-Dyson random matrix phenomena. 
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This principle will be explained in the present section. We mention in passing 
that a fascinating recent development (Kit08, SRF09] uses the same principle 
for a homotopy classification of topological insulators and superconductors. 



3.1 Fock space setting 

We now describe an extended setting, which replaces the Hermitian vector 
space V by its exterior algebra A(V) but otherwise retains Dyson's setting to 
the fullest extent possible. In physics language we say that we pass from the 
(single-particle) Hilbert space V to the fermionic Fock space A(V) generated 
by V. The Z-grading A(V) = ©„ A™ (V) by the degree n has the physical 
meaning of particle number. Thus A°(V) = C is the vacuum, A-^V) = V is the 
one-particle space, f\ 2 {V) is the two-particle space, and so on. We adhere to 
the assumption of finite-dimensional V . Particle number n then is in the range 

< n < N := dim V. Note that dim A n (V) 



n 

The n-particle subspace A n (V) of the Fock space of a Hermitian vector space 
V carries an induced Hermitian scalar product defined by 



(ui A • • • A u n , vi A • • • A «n) A «(v) = Dct 



({ui,vi) v ... {ui,v n ) v ^ 



\{u n ,Vl)v ■■■ {U 



• (3-1) 



Relevant operations on A(V) are exterior multiplication (or particle creation) 
e(v) : A n (V) — > A n+1 (V) by v G V and contraction (or particle annihilation) 
t(f) : A n (V) — > A n ~ 1 {V) by / G V*. The standard physics convention is to 
fix some orthonormal basis {ek}k=i,...,N of V and write at := e(efc) for the 
particle creation operators and := i((efc,-)y) for the particle annihilation 
operators. This notation reflects the fact that Hermitian conjugation f in Fock 
space relates e(v) and t(f) by e(v)^ = t(f) where / = (v, -)y. The operators a\ 
and afc satisfy the so-called canonical anti-commutation relations 

a k ai + aia k = = a\a\ + a\a\ , a[a l + a^l = S M . (3.2) 

These represent the defining relations of the Clifford algebra C1(W) of the vector 
space W = V © V* with quadratic form (v ©/,«'© /') h-> f(v') + f'(v). 

Having introduced the basic Fock space structure, we now turn to what is 
going to be our definition of a symmetry group G acting on Fock space A(V). 
As before, we assume that we are given a normal subgroup Go C G. The action 
of the elements g G Go is defined by unitary operators on V which are extended 
to A(V) by 

g(vi A • • • A v n ) := {gvi) A • • • A (gv n ) . (3.3) 

Similarly, the anti-unitary operator of time reversal T is defined on V and is 
extended to A(V) by T(vi A ■ ■ ■ A v n ) := (Tv x ) A ■ ■ ■ A {Tv n ). 
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Now the Z-grading of Fock space offers the natural option of introducing 
another kind of anti-unitary operator, which is a close cousin of the Hodge star 
operator for the deRham complex: particle-hole conjugation, G, transforms an 
n-particle state into an (N — n)-particle state or a state with n holes. (Note 
the change of meaning of the symbol G as compared to Section [3.31 ) 

Definition 3.1 Fix a generator Ct £ A N (V), N = dim^, with normalization 
(Q,Q) A Nfv) = 1- Then particle-hole conjugation C : A n (V) — > A N ~ n (V) is the 
anti-unitary operator defined by 

(C^)A^ = (^'>a»(V)«- 

Thus the definition of the operator G uses the Hermitian scalar product of Fock 
space and a choice of fully occupied state Cl. An elementary calculation shows 
that C 2 \ A n (v) = (-l) n ( N ~ n ). 

What are the commutation relations of G with T and g £ Gq ? To answer 
this question, we observe that by dim (V) = 1 we may always choose Ct to 
be T-invariant (i.e., TCI = Cl). Then from the following computation, 

(rcvo a Tip' = t ((ap) a ip') = t ((v>, v') n) 

= {ip,ip')TCl = (T^,Tip')n = (CTip) A Tip', 

we have CT = TC. Also, making the natural assumption that both the vacuum 
space A°(V) and the fully occupied space A N (V) transform as the trivial Go- 
representation (i.e., g£l = for g G Go), a similar calculation gives Cg = gC. 

In order to enlarge the set of possible symmetries and hence the scope of 
the theory, we now introduce a 'twisted' variant of particle-hole conjugation. 
Let S £ U(F) be an involution (S 2 = Id) and extend S to A(V) by Eq. (HP) . 
To obtain an extension of the group Go , we require that S commutes with T, 
satisfies SCI = CI, and normalizes Go , i.e., SGqS -1 = Go • (Here we identify Go 
with its action on Fock space.) By twisted particle-hole conjugation we then 
mean the operator G = CS = SC. Note that GG G _1 = G and CT = TC. 

Definition 3.2 On the Fock space A(V) over a Hermitian vector space V , let 
there be the action of a group G = Go U TGo U GGo U CTGq with Go a normal 
subgroup and CT = TC. We call this a minimal extension of Dyson's setting 
if Go acts by unitary operators defined on V , time reversal T acts as an anti- 
unitary operator also defined on V , and twisted particle-hole conjugation C is 
an anti-unitary bijection A n (V) — > A N ~ n (V). 

3.2 Classification goal 

The simplest question to ask now is this: what is the structure of the set of 
Hamiltonians that operate on Fock space A(V) and commute with the given 
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G-action on A(V)? Since this question ignores the grading of Fock space by 
particle number, it takes us back to Dyson's setting and the answer is, in fact, 
provided by Dyson's Threefold Way. (Note that in the absence of restrictions, 
the most general Hamiltonian in Fock space involves n-body interactions of 
arbitrary rank n = 1, 2, 3, . . . , N.) So there is nothing new to discover here. 

We shall, however, be guided to a new and interesting answer by asking a 
somewhat different question: what is the structure of the set of one-body time 
evolutions of A(V) which commute with the given G-action? Here by a one- 
body time evolution we mean any unitary operator obtained by exponentiating 
a self- adjoint Hamiltonian H which is quadratic in the particle creation and 
annihilation operators: 

H = J2 ki Wkl 4°i + g J2ki ( Zkl a k a l + Zkl aja fc ) . (3.4) 

These operators U = e~ ltH / h g U(AV) form what is called the spin group, 
Spin(Wk), of the 2iV-dimensional Euclidean R- vector space Wr spanned by the 
Majorana operators a k + a' k , \a k — iat (k = 1, . . . , N). Spin(Vt / ig) ~ Spin 2 7v is a 
double covering of the real orthogonal group SO(Wr) ~ SO2N . The spin group 
of most prominence in physics is Spin 3 = SU2 , a double covering of SO3 . (This 
double covering is known to physicists by the statement that a rotation by 2tt, 
which acts as the neutral element of SO3 , changes the sign of a spinor.) 
Thus, our interest is now in the set 

Z Spin (G) := Z V{AV) (G) n Spm(W R ) (3.5) 

of G-compatible time evolutions in Spin(VF]g) C U(AV). By adaptation of the 
earlier definition ()2.4f) . the G-compatibility conditions are 

U = g^Ug^ 1 = gi U- l g^ (3.6) 

for all g G G U CTG and gi G TG U CG . 

3.3 Reduction to Nambu space 

To investigate the set Zs p i n (G) we use the following fact. Any invertible element 
A G C1(W) determines an automorphism 7 1— > A'yA' 1 of the Clifford algebra 
C1(W) by conjugation. This conjugation action restricts to a representation 

r(g)w := gwg' 1 (3.7) 

of Spin(WR) C Cl(VF) on Wj& C C1(W). Phrased in physics language, the 
set of Majorana field operators w = J2k( z k a k + ^fc a fc) ^ i s closed under 
conjugation w \-> gwg^ 1 by one-body time evolution operators g G Spin(WR). 
In fact, by elementary considerations one finds that r(g) : w 1— > gwg -1 for 
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g G Spin(W / K) is an orthogonal transformation r(g) G SO(W / ir) of the Euclidean 
vector space Wr ■ The mapping r : Spin(WR) — > SO(Wr), g i-> r(p) = r(—g) 
is two-to-one. It is a covering map, which amounts to saying that any path in 
SO(Wr) lifts uniquely to a path in SpinfWk). Note also that the linear mapping 
r(g) : Wr — > Wr extends to a linear mapping -r(g) : W — > W by C-linearity. 

Thus, instead of studying Spin(WiO as a group of operators on the full 
Fock space A(V), we may simplify our work by studying its representation 
r : Spin — > SO on the smaller space W = V © V* , here referred to as Nambu 
space. Of course the object of interest is not Spin(WR) but its intersection with 
the G-compatibility conditions. To keep track of the latter conditions, we now 
transfer the G-action from A(V) to W = V © V*. 

It is immediately clear how to transfer the actions of Go and T, as these are 
defined on V . In the case of the twisted particle-hole conjugation operator G, 
we do the following computation. Let rp G A n (V) and ip' G A n+ (V). Then 

(g4v) a = (Sa{i;, v>') n = (Sip, SakS- 1 ^') n 

= {CiP) A (SakS- 1 )^' = {-^-^(SakS^CiP) A t// . 

Thus the twisted particle-hole conjug ate ofa|eFc Cl(W) is Ca\C~ x = 
±5'ofc5' _1 G V* C C1(W) where the sign alternates with particle number. Note 
that the operation at i— > ±ScikS~ 1 is anti-unitary. Note also that the untwisted 
particle-hole conjugation at i— > a k is none other than the C-antilinear bijection 
Cy : V^V*,v^ (v,-) v . 

To sum up, we have the following induced structures on Nambu space: 

• One-body time evolutions g G Spin(WR) act on W = V@V* by orthogonal 
transformations r(g) G SO(Wr). 

• G is defined on V and acts on W = V © V* by g(v © /) = gv © {g~ 1 ) 1 f '. 
The same goes for time reversal T. 

• The operator C of twisted particle- hole conjugation induces on W = 
V © V * an anti- unitary involution V V* . 

The goal of symmetry classification now is to characterize the set Zso(G) of 
elements in SO(Wr) which are compatible with the induced G-action on W. 
This problem was posed and solved in [HHZ05] . by using an elaboration of the 
algebraic tools of Section [2] to make a reduction to the case of the trivial group 
Go = {Id}. (The involution V V* given by twisted particle- hole conjugation 
is called a mixing symmetry in [HHZ05J.) The outcome is as follows. 

Theorem 3.1 The space Zso(G) is a direct product of factors each of which 
is a classical irreducible compact symmetric space. Conversely, every classical 
irreducible compact symmetric space occurs in this setting. 
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There is no space to reproduce the proof here, but in order to turn the theorem 
into an intelligible statement we now record a few basic facts from the theory 
of symmetric spaces |Hel78|, ICM04] . 



3.4 Symmetric spaces 

Let M be a connected m-dimensional Riemannian manifold and p a point of 
M. In some open subset N p of a neighborhood of p there exists a map s p : 
N p — > N p , the geodesic inversion with respect to p , which sends a point x G N p 
with normal coordinates [x\ , . . . , x m ) to the point with normal coordinates 
(— xi, . . . , —x m ). The Riemannian manifold M is called locally symmetric if 
the geodesic inversion is an isometry (i.e., is distance-preserving), and is called 
globally symmetric if s p extends to an isometry s p : M — > M, for all p G M. A 
globally symmetric Riemannian manifold is called a symmetric space for short. 

The Riemann curvature tensor of a symmetric space is covariantly constant, 
which leads to three distinct cases: the scalar curvature can be positive, zero, 
or negative, and the symmetric space is said to be of compact type, Euclidean 
type, or non-compact type, respectively. In random matrix theory each type 
plays a role: the first one provides us with the scattering matrices and time 
evolutions, the second one with the Hamiltonians, and the third one with the 
transfer matrices. Our focus here will be on compact type, as it is this type 
that houses the unitary time evolution operators of quantum mechanics. 

Symmetric spaces of compact type arise in the following way. Let U be 
a connected compact Lie group equipped with a Cartan involution, i.e., an 
automorphism r : U — > U with the involutive property r 2 = Id. Let K C U 
be the subgroup of r-fixed points u = t(u). Then the coset space U/K is a 
compact symmetric space in a geometry defined as follows. Writing u := Lie(£/) 
and 6 := Lie(-fT) for the Lie algebras of the Lie groups involved, let u = I © p 
be the decomposition into positive and negative eigenspaces of the involution 
dr : u — > u induced by linearization of r at unity. Fix on p a Euclidean scalar 
product (-,-) p which is invariant under the adjoint i'T-action Ad(k) : p — > p 
by I 4 kXk^ 1 . Then the Riemannian metric g u x evaluated on vectors v,v' 
tangent to the coset uK is g u K{v, v') := {dL~ l {v) : dL~ 1 (u'))p where dL u denotes 
the differential of the operation of left translation on U/K by u G U. 

It is important for us that the coset space U/K can be realized as a subset 



by the Cartan embedding U/K M C U, uK \— > ut(u~ 1 ). The metric tensor 
in this realization is given (in a self-explanatory notation) by g = Tr dx dx^ 1 . It 
is invariant under the X-action M — > M by twisted conjugation x i— > uxt^u^ 1 ). 
The geodesic inversion with respect to y G M is s y : M — > M, x t-t yx^ 1 y. 




(3.8) 
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family 


compact type 


Euclidean type 




A 




H complex Hermitian 




Al 




H real symmetric 




All 


U 2 Ar/USp 2iv 


H quaternion self-dual 




C 


USp 27V 


Z complex symmetric 




CI 


USp 2Ar /UAr 


Z complex sym., W = 





B,D 




Z complex skew 




Dill 


S0 2 tv/Uat 


Z complex skew, W = 





Alll 


IW(U P x u g ) 


Z complex p x q, W = 





BDl 


p+q /(0 p x 0,) 


Z real p x q, W = 




Cll 


USp 2p+2g /(USp 2p x USp 2(? ) 


Z quaternion 2p x 2g, 


W = 



Table 1: The Cartan table of classical symmetric spaces 



We note that special examples of compact symmetric spaces are afforded 
by compact Lie groups K. For these examples, one takes U = K x K with flip 
involution r(k, k') = (A/, k) leading to U/K = (K xK)/K ~ K. Cartan's list of 
classical (or large families of) compact symmetric spaces is presented in Table HJ 
The form of the generator H of time evolutions indicated in the 

third column, where the notation W, Z refers to the Fock space expression (|3.4p 

[W Z \ 

which translates to H = I ^ ^ t j by the covering map r : Spin — > SO . 
3.5 Post-Dyson classes 

We now run through the symmetry classes beyond those of Wigner-Dyson. As 
was mentioned before, these appear in various areas of physics and in the ran- 
dom matrix theory of L-functions. For brevity we concentrate on their physical 
realization by quasi-particles in disordered metals and superconductors. 

3.5.1 Class D 

Consider a superconductor with no symmetries in its quasi-particle dynamics, so 
G = {Id}. (Some concrete physical examples follow below.) The time evolutions 
u = e~ ltH / h in this case are constrained only by the condition u £ Spin(Wig) 
in Fock space and r(u) E SO(Wr) in Nambu space. The orthogonal group 
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SO(VKr) ~ S02AT is a symmetric space of the D family - hence the name class 
D. In a basis of Majorana fermions a k + a, , \a k — iat , the matrix of \H G SO2N 
is real skew, and that of H is imaginary skew. 

Concrete realizations are found in superconductors where the order param- 
eter transforms under rotations as a spin triplet in spin space and as a p-wave 
in real space. A recent candidate for a quasi-2d (or layered) spin-triplet p- 
wave superconductor is the compound S^Ru O4 . (A non-charged analog is 
the A- phase of superfluid 3 He.) Time-reversal symmetry in such a system may 
be broken spontaneously, or else can be broken by an external magnetic field 
creating vortices in the superconductor. 

The simplest random matrix model for class D, the SO-invariant Gaussian 
ensemble of imaginary skew matrices, is analyzed in Mehta's book |Meh04| . 
From the expressions given there it is seen that the level correlation functions 
at high energy coincide with those of the Wigner-Dyson universality class of 
unitary symmetry (Class A). The level correlations at low energy, however, 
show different behavior defining a separate universality class. This universal 
behavior at low energies has immediate physical relevance, as it is precisely the 
low-energy quasi-particles that determine the thermal transport properties of 
the superconductor at low temperatures. 

3.5.2 Class Dili 

Now, let time reversal T be a symmetry: G = {Id, T}. Physically speaking this 
implies the absence of magnetic fields, magnetic impurities and other agents 
which distinguish between the forward and backward directions of time. Our 
physical degrees of freedom are spin- 1/2 particles, so T 2 = —Idy ■ 

According to (I3.6P we are looking for the intersection Zso(G) of the condi- 
tion = TuT" 1 with Spin(W]R), or after transfer to Nambu space, SO(Wr). 
By introducing the involution t(u) := TuT~ l we express the wanted set as 

^so(G) = {u G SO(Wr) I u- 1 = t(u)} . (3.9) 

Following the discussion around Eq. (|3.8p we have Zso(G) ~U/K where U = 
SO(Wik) and K C U is the subgroup of r- fixed points in U. 

In order to identify K we note that time reversal T : W —> W preserves the 
real subspace of Majorana operators a. + at , ia fc — \a\ . Because T 2 = —Id, 
the R-linear operator T : Wr — > W$t is a complex structure of the real vector 
space ~ M? N ~ C N . In other words, there exists a basis {eij , e2,j}j=i,...,N 
of Wr such that Teij = e2j and T&i,j = ~ e i,j ■ Now the r-fixed point 
condition k = r(k) says that k G K commutes with the complex linear extension 
J : W — > W of T : Wr — > Wr by Je±j = ±ie±,j where e±j = eij ± ie2j ■ 
The general element k with this property is a U^-transformation which acts 
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on span c {e +i i, . . . , e +j at} as k and on span c {e_ ; i, . . . ,e_ ; ;v} as k = (k 1 ) t . 
Hence K = Ujv and 



which is a symmetric space in the Dili family. 

Known realizations of this symmetry class exist in gapless superconductors, 
say with spin-singlet pairing, but with a sufficient concentration of spin-orbit 
impurities to break spin-rotation symmetry. In order for quasi-particle excita- 
tions to exist at low energy, the spatial symmetry of the order parameter should 
be different from s-wave. A non-charged realization occurs in the S-phase of 
3 He, where the order parameter is spin-triplet without breaking time-reversal 
symmetry. Another candidate are heavy-fermion superconductors, where spin- 
orbit scattering often happens to be strong owing to the presence of elements 
with large atomic weights such as uranium and cerium. 

3.5.3 Class C 

Next let the spin of the quasi-particles be conserved, but let time-reversal sym- 
metry be broken instead. Thus magnetic fields (or some equivalent T-breaking 
agent) are now present, while the effect of spin-orbit scattering is absent. The 
symmetry group of the physical system then is G = G$ = Spin 3 = SU2 . Such a 
situation is realized in spin-singlet superconductors in the vortex phase. Promi- 
nent examples are the cuprate superconductors, which are layered and exhibit 
an order parameter with d-wave symmetry in their copper-oxide planes. 

The symmetry-compatible time evolutions are identified by going through 
the process of eliminating the unitary symmetries Go = G. For that, we de- 
compose the Hilbert space asV = E®R,E = HomdR, V), where R := C 2 
is the fundamental representation of G = SU2 . Now there exists a skew- 
symmetric SU2-equivariant isomorphism (known in physics by the name of 
spin-singlet pairing) between the vector space R and its dual R* . Therefore 
we have W = V @V* ~ (E © E*) © R , and elimination of the conserved factor 
R transfers the canonical symmetric form of W = V © V* to the canonical 
alternating form (e © /, e' © /') i-> f(e') - /'(e) of E © E*. On transferring 
also the Hermitian scalar product from V © V* to E © E* , one sees that the 
G-compatible time evolutions form a unitary symplectic group, 



which is a compact symmetric space of the C family. 
3.5.4 Class CI 

The next class is obtained by taking spin rotations as well as the time reversal 
T to be symmetries of the quasi-particle system. Thus the symmetry group 



Zso(G) ~ U/K ~ SCWUtv , 



(3.10) 



Z so (G) = SO(W R ) G ~ USp(£ © E*) , 



(3.11) 
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now is G = Go U TGq with Go = Spin 3 = SU2 . Like in the previous symmetry 
class, physical realizations are provided by the low-energy quasi-particles of 
unconventional spin-singlet superconductors. The superconductor must now be 
in the Meissner phase where magnetic fields are expelled by screening currents. 

To identify the relevant symmetric space, we again transfer from V © V* = 
(E © E*) <S> R to the reduced space E © E*. By this reduction, the canonical 
form undergoes a change of type from symmetric to alternating as before. We 
must also transfer time reversal; because the fundamental representation R = 
C 2 of SU2 is self-dual by a skew-symmetric isomorphism, the parity of the 
time-reversal operator changes from T 2 = — Idy^v to T 2 = +Id£©.E* by the 
mechanism explained at the end of Section 12.21 

We have Z so (G) ~ U/K where U = USp(£ , ffi£'*) and K is the subgroup of 
elements fixed by conjugation with T. Because the reduced T squares to +1, we 
may realize it on matrices as the complex conjugation operator by working in a 
basis of T-fixed vectors of E © E* . The Lie algebra elements X G usp(E © E*) 



have the form X = I — —r I with anti-Hermitian A and complex symmetric 



B. They commute with the operation of complex conjugation if A is real skew 
and B real symmetric. Matrices X with such A and B span the Lie algebra 
Ujv, N = dim(E). At the Lie group level it follows that K = \] n . Hence 
Zso(G) ~ USp2Ar/Ujv - a symmetric space in the CI family. 

3.5.5 Class AIII 

So far, we have made no use of twisted particle-hole conjugation C as a sym- 
metry, but now let the symmetry group be G = Go U CGo where Go = Ui acts 
on W = V © V* by v © / ^ zv © z~ l f (for z G C , \z\ = 1). 

In order for the elements u £ SO(Wr) to commute with the Go-action, 
they must be of the block-diagonal form u = k © (/c -1 )', k € U(V). Therefore 
^so(Go) ^ U(V). The wanted set then is Z so (G) ~ U/K with U = U(V) and 
K the subgroup of elements which are fixed by conjugation with C. 

Recall from Section ETT1 that C\y = CS where S G U, S 2 = Id, and un- 
twisted particle-hole conjugation G coincides (up to an irrelevant sign) with 
the canonical bijection Cy : V — > V*, Cy(v) = (v,-)v- The condition for 
u = k © {k~ 1 ) 1 to belong to K reads 



Since k" 1 = k' and C~ 1 k t C = k', this condition is equivalent to k = SkS . 

Now let V = V+ © V- where V± are orthogonal subspaces with projection 
operators II-|- . Then if S = II + — II_ we have K = U(V+) x U(VL) and hence 




(k-y = ckc- 1 . 



(3.12) 



Z so (G) ~ U(V)/(\J(V + ) x \J(y_)) 



(3.13) 
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or Zso(G) — Uat/(U p x Uat_ p ) with p = dimy + . 

The space (|3.13p is a symmetric space of the j4III family. Its symmetry class 
is commonly associated with random-matrix models for the low-energy Dirac 
spectrum of quantum chromodynamics with massless quarks [Ver9 4] . An al- 
ternative realization exists [ASZ02J in T-invariant spin-singlet superconductors 
with d-wave pairing and soft impurity scattering. 

3.5.6 Classes BDI and CII 

Finally, let the symmetry group G have the full form of Definition 13.21 with 
Co = Ui and C as before (Class AQI) and a time-reversal symmetry T, T 2 = 
±Id. We recall that the elements of Z S o(G ) are u = k © k £ U(V). 

The requirement of commutation with the product (j> := CT : V — > V* of 
anti-unitary symmetries is equivalent to the condition <\> = k l $k. 

Let U := Z S o(G UCTG ). To identify U, we use that <j)(v)(v') = {STv,v') 
and ST = TS. By the computation of (|2.9|) . it follows that the parity of T 
equals the parity of the isomorphism <f> : V — > V*. In other words, if T 2 = eld 
then (fr = ecf). Thus the condition (ft = k <p k singles out an orthogonal group 
U = O(V) — On in the symmetric case (e = +1) and a unitary symplectic 
group U = USp(y) ~ USp^y in the alternating case (e = — 1). 

In both cases, the wanted set is Zso (G) = U/K with K the subgroup of fixed 
points k = C' _1 (fc _1 )'C' = SkS. In the former case we have K ~ O p x On- p , 
and in the latter case K ~ USp p x USpjv_ p (with even N, p). Thus we arrive 
at the final two entries of Cartan's list: 



These occur as symmetry classes in the context of the massless Dirac operator 
|Ver94] . Class BDI is realized by taking the gauge group to be SU2 or USp 2 „ , 
Class CII by taking fermions in the adjoint representation or gauge group SO n . 

4 Discussion 

Given the classification scheme for disordered fermions, it is natural to ask 
whether an analogous scheme can be developed for the case of bosons. Although 
there exists no published account of it (see, however, [LSZ 06]). we now briefly 
outline the answer to this question. 

The mathematical model for the bosonic Fock space is a symmetric algebra 
S(V). It is still equipped with a canonical Hermitian structure induced by 
that of V. The real form of Nambu space W = V © V* for bosons has an 
interpretation as a classical phase space spanned by positions qj = (a - +0.^)7^/2 



Class BDI : U/K ~ N /(0 P x N _ p ) (T 2 
Class CII : U/K ~ US Pjv /(USp x \JSp N ) (T 2 



+1), 
-!)• 




At the level of one-body unitary time 
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evolutions in Fock space, the role of the spin group Spin(WjR) for fermions is 
handed over to the metaplectic group Mp(Wr) for bosons. 

By the quantum-classical correspondence, a one-parameter group of time 
evolutions ut = e~ ltH / h g Mp(Wr) in Fock space gets assigned to a linear 
symplectic flow r(u t ) G Sp(Wr) in classical phase space. This correspondence 
r : Mp(Wr) — > Sp(W / k) is still two-to-one (reflecting, e.g., the well-known 
fact that the sign of the harmonic oscillator wave function is reversed by time 
evolution over one period). An important difference as compared to fermions 
is that the classical flow r(u t ) G Sp(Wk) is not unitary in any natural sense. 

In nuclear physics, the differential equation of the flow r(ut) is called the 
RPA equation. For example, in the case without symmetries this equation reads 

J t 4 = J2j ( a W + a i B i k ) ' J t ak = ^ jCjk + ajDjk ^ ' ^ 

where one requires B = B 1 , C = C*, and D = —A 1 in order for the canonical 
commutation relations of the boson operators a^,a to be conserved. Unitarity 
of the flow (as a time evolution in Fock space) requires A = —A^ and C = B^ . 
This should be compared with the fermion problem in Class C, where one has 
exactly the same set of equations but for a single sign change: C = —B*. Thus 

A B^ 



the corresponding generator of time evolution is X = ^_|_^> ^[J where the 

plus sign applies to bosons and the minus sign to fermions. In either case X 
belongs to the same complex Lie algebra, sp(PF). The difference is that the 
generator for fermions lies in a compact real form usp(I^) C sp(VF), whereas 
the generator for bosons lies in a non-compact real form sp(Wk) C sp(PF). 

This remains true in the general case with symmetries. Thus if the word 
'symmetry class' is understood in the complex sense, then the bosonic setting 
does not lead to any new symmetry classes; it just leads to different real forms 
of the known symmetry classes viewed as complex spaces. The same statement 
applies to the non-Hermitian situation. Indeed, all of the spaces of [Mag08 
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complex or non-compact real forms of the symmetric spaces of Cartan's table. 
Here we must reiterate that the notion of symmetry class is an algebraic one 
whose prime purpose is to inject an organizational principle into the multitude 
of possibilities. It must not be misunderstood as a cheap vehicle to produce 
immediate predictions of eigenvalue distributions and universal behavior! 

Let us end with a few historical remarks. The disordered harmonic chain, 
a model in the post-Dyson Class BDI, was first studied by Dyson |Dys53| . 
The systematic field-theoretic study of models with sublattice symmetry (later 
recognized as members of the chiral classes AQI, BDI) was initiated by Op- 
permann, Wegner and Gade |OW79t IGad931 ICW91j . Gapless superconductors 
were the subject of numerous papers by Oppermann; e.g., [Opp90| computes 
the one-loop beta function of the non-linear sigma model for Class CI. 
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The 10- way classification of Section [3] was originally discovered by a very 
different reasoning: the mapping of random matrix problems to effective field- 
theory models [Zir96j combined with the fact that closure of the renormalization 
group flow takes place for non-linear sigma models where the target is a sym- 
metric space. A less technical early confirmation of the 10-way classification 
came from Wegner's flow equations [Weg94~] . These take the form of a double- 
commutator flow for Hamiltonians H belonging to a matrix space p; if the 
double commutator [p, [p,p]] closes in p, so does Wegner's flow. The closure 
condition is satisfied precisely if p is the odd part of a Lie algebra u = t © p 
with involution, i.e., the infinitesimal model of a symmetric space. 

Last but not least, let us mention the viewpoint of Volovik (see, e.g., fVol03]) 
who advocates classifying single-particle Green's functions rather than Hamil- 
tonians. That viewpoint in fact has the advantage that it is not tied to non- 
interacting systems but offers a natural framework in which to include (weak) 
interactions. 
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